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Abstract: 

The semileptonic decays of polarised top quarks are analysed for a general chirality- 
conserving tbW vertex. We calculate double differential distributions for the charged 
lepton and the neutrino to order a s in the QCD coupling. We present these QCD 
corrections in terms of compact parameterisations that should be useful for the 
future investigation of the structure of the top decay vertex on the basis of large 
data samples. 
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So far the experimental information about the decays of top quarks is not very 
detailed. Data from the Tevatron are consistent with the expectation that these 
decays are governed by the V — A charged current interactions of the standard 
model (SM); but because of the size of the present experimental errors sizable new- 
physics effects can not be excluded (c.f. for instance [I]). 

In view of the extremely large mass of the top quarks and the circumstance 
that they do not hadronise these particles are excellent probes of new interactions 
that may uncover at energies of a few hundred GeV. There is a vast literature on 
the phenomenology of such interactions in top quark decays (for overviews, see for 
instance IS])- One possibility is a small V + A admixture to the standard left- 
handed t — > b current as is predicted for instance by SU{2) L x SU{2) R x U(l) 
extensions of the standard model. The measured branching ratio of the radiative 
weak decay b — ► 37 provides a stringent constraint jU El El HI on the coupling strength 
Kr of such a right-handed admixture: \kr\ < 0.04. Of course, one may envisage 
contrived scenarios where other contributions to b — > 37 lead to cancellations that 
invalidate this bound - in any case, a direct search for such a coupling in t and t 
decays is indispensable. Sensitive tests require polarised top quarks. At the LHC 
the single top production processes will yield large samples of polarised t and i 
quarks, and it has been estimated [HI E] that a statistical sensitivity 5kr ~ 0.06 
can be reached. From tt pair production at a high-luminosity linear e + e~ collider a 
sensitivity 5kr ~ 0.03 will be feasible [TU1 ITU 13]. 

Energy and angular distributions for polarised top quark decays with a V + A 
admixture in the tbW vertex were investigated in a number of articles, including 
O H31 UJJ [JS] . In ref . [T3| semileptonic t decays were analysed and it was pointed 
out that the polarisation-dependent part of the neutrino distribution (which can be 
determined experimentally by measuring the missing momentum) is very sensitive to 
Kr. The order a s QCD corrections for the V — A part of the charged current [161 177] 
were incorporated in [TjJ]. In view of the above-mentioned expected sensitivity to 
hr at future colliders one should note that taking into account QCD corrections in 
future data analysis is mandatory, because they can mimic a small V + A admixture. 

In this paper we extend the work of J3] i n that we compute the order a s QCD 
corrections to both the V — A and V + A Born contributions to t — > biv?. Moreover 
we take the finite width of the intermediate W boson and the non-zero mass of the b 
quark into account. The latter is important in the case of a small V + A admixture, 
because terms linear in kr in the lepton distributions, which we compute below, 
require a chirality flip of the b quark. For small kr, say \k,r\ < 0.1, the terms 
proportional to URin^/mt are of the same order of magnitude as the contributions 
proportional to k r . Furthermore we compute a T-odd triple correlation which is 
generated by a non-standard CP- violating phase that may be present in the general 
chirality-conserving tbW vertex. 

We consider the semileptonic decay of spin-polarised top quarks to order a s in 
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the QCD coupling, which amounts to studying the reactions 



t -> Hp&)+^ + (p*)+^) + <?(p s )- (1) 

All momenta refer to the rest frame of the top quark. As far as the tbW vertex is 
concerned we use the generalised interaction 

C tb = -^=V tb iY(a-[3 l5 )bW ll + h.c, (2) 

where a, (5 are complex couplings. A right-handed admixture to the V — A charged 
current interaction of the SM (a, — (5 — 1) will be parameterised by choosing 
a = 1 + kji, (3 — 1 — Kr. In the calculations below we neglect the lepton masses but 
take the mass of the b quark and the finite width of the intermediate W boson into 
account. 

In the calculation of the differential decay distributions for ([TJ we have performed 
the quark wave function renormalisations that remove the ultraviolet divergencies 
in the on-shell scheme. The infrared divergencies are canceled using a standard 
phase space slicing procedure. The phase space of the four-particle final state in the 
reaction t biv^g is split into two disjoint regions where the (scaled) energy x g = 
2E g /m t of the gluon is smaller and larger than an arbitrary, but small separation 
parameter x m i n , respectively: 1 = 6(x m j n — x g ) + Q(x g — x m i n ). Integrating the 
respective squared matrix element over the phase space of the soft gluon (x g < x m i n ) 
and adding the result to the order a s squared matrix element for the three particle 
final state b£vg yields the infrared finite differential decay distribution dT s+v+soft 
(B = Born, V = virtual) and dThard, which describes radiation of "resolved" gluons 

{Xg J> X m i n ^j. 

For checks of possible deviations from the V — A law in the tbW vertex useful 
observables are the double differential energy-angle distributions dT/dxid cos 9i (i = 
£,u) for the charged lepton and for the neutrino, where Xi = 2Ei/m t and O-i is 
the angle between the three-momentum of the lepton % and the unit vector s that 
specifies the polarisation direction of the ensemble of top quarks in the t rest frame. 
The degree of polarisation is denoted by S = |s|. Adding the contributions from 
dT B+v+soft and from dT hard , these distributions can be put into the following form: 



dT it/ ,20 m t 

" \ V tb\ — - — 



dxtd cos 6 



\Vtb\ o 



doc it d cosy, 



(F* n + C F — Ff n ) + (F* s + C F — Ff s )S cos9 t , (3) 

7T 7T 



(F™ + C F — F» n ) + {F» S + C F — F? S )S cosOv , (4) 

7T 7T 



where C F = 4/3 and < Xg,x u < (mf — ml)/mf. The corresponding distributions 
for the anti-top quark are obtained by changing the sign of the term proportional 
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to cos 6{. From these distributions the corresponding 2x2 top decay spin density 
matrices can be extracted in straightforward fashion 1 . 

The (i = £, is, j = n, s, and a = 0, 1) which are functions of xe and x u , 
respectively, can be decomposed according to the contributions from the vector and 
axial vector couplings a, (3 : 

K 3 = W?Ft + Re (a*(3)F^ + \(3\ 2 F^ . (5) 

The following relations hold between the functions that appear in the charged lepton 
distribution and those of the corresponding neutrino distribution: 

r a,a a,a ' a,P a,P ' a,a(3 r a,a(3 ' 

rpls jpi/s Tpis -pus Tpls jpus (p\ 

V" r a,ai r a,f3— r a,/3 ) r a,a/3 ~ r a,a/3 > \ v J 

where a = 0, 1. For the sake of presenting compact expressions we write down the 
lowest order functions of the charged lepton distribution in terms of one-dimensional 
integrals (r = a,/3,a/3): 

F °' r= 512^ dx b D w (x b )^ r (x e ,x b ), (7) 

where D w = ((1 - x b - m 2 w + mf) 2 + m^f 2 ^) -1 , m w = m w /m t , T w = T w /m t , 
m b = m b /m t , x^ min = (1 - x e + m b /(l - xi)), x b , max = (1 +m b ), and 

Ffc = -2m 3 b +m 2 (x b -2) + 2m b (x b -l)-x 2 b + x b (3-2x £ )-2(x e -l) 2 , 

Ffctf = 2(1 + m\ - x b ){2 - x b - 2x e ) , 

Ffy = 2ml + m 2 {x b -2)-2m b (x b -l)-x 2 b +x b (3-2x £ )-2(x e -l) 2 1 

F o% = -(! + m l ~ x b )(2m 2 b + x b (x e - 2) + 2(x e - l) 2 + 2m b x e ) , 

F l Q * af} = 4(1 + m\ - x b ) 2 + 2(1 - m\ - x b ) (x b - A)x e - A(x h - 2)x\ - 4rrf , 

F o,p = H± + ™l - Xb)(2rh 2 b + x b (x e - 2) + 2(x e - l) 2 - 2m b x e ) . (8) 

The corresponding functions for the neutrino distribution are obtained with the 
relations (jUJ). Using a, (3 = 1 ± k,r it is easy to see that those terms in the Born 
distributions which are linear in k,r are accompanied by a factor m b (or ml). 

The functions F[ 3 a , Fl 3 af3 , and F^a which are induced by the order a s corrections 
are shown in Figs. 1 and 2. We use m t = 174.3 GeV, m b = 4.8 GeV, mw = 80.42 
GeV, and Tw = 2.12 GeV. For SM couplings a — (3 — 1 and for the narrow width 
approximation Tw — > we can numerically compare our results with those of 
and we find agreement. 

The size of the QCD corrections is typically 6 to 8 percent for Xi < 0.9 and 
somewhat larger at the upper end of the spectrum. For Xi — > 1 the curves show the 

1 These density matrices enter the calculations of the (differential) cross sections that describe 
t and/or i production and decay in the on-shell approximation. 



3 

















A [10- 7 ] 


-8.06 


-10.54 


-14.37 


-7.55 


-9.86 


-15.25 




5.50 


5.50 


4.96 


5.24 


5.24 


5.18 


Bo 


-0.2853 


-0.3087 


-1.6414 


-1.1058 


-1.2267 


-1.6760 


Bi 


3.3781 


3.6545 


22.6920 


17.7909 


19.7949 


23.4286 


B 2 


-16.3481 


-17.902 


-131.1934 


-122.4657 


-136.8286 


-136.6096 


B 3 


40.0362 


44.3232 


408.7910 


467.4425 


523.4506 


428.1320 


B 4 


-53.1628 


-59.4165 


-741.6436 


-1083.4018 


-1213.9973 


-781.1829 


B 5 


36.4636 


41.0886 


784.5948 


1563.7047 


1751.0343 


831.2825 


B 6 


-10.1118 


-11.4767 


-448.9931 


-1374.8811 


-1536.9986 


-478.5832 


B 7 






107.4585 


674.4979 


752.1898 


115.2665 


B 8 








-141.6121 


-157.4561 





Table 1: The coefficients which determine the fits (JUJ) to the order a s QCD contri- 
butions to the charged lepton distribution (jSJ). A an d A\ are given in units of 10~ 7 
and 10~ 5 , respectively. 

developoment of the well-known logarithmic singulariy [16 which can be removed 
by exponentiation. Rather than presenting the exact analytical formulae for the 
functions F^ a , F[ 3 a3) and F^g, which are quite lengthy, we give them in terms of 
simple parameterisations determined by fits. We required the fits to be better than 
4% in the range < xg, x u < 0.98. This requirement can be satisfied by using 
a rational function for xi,x u < 0.218 and a polynomial of at most 8th order for 
xe,x u > 0.218. This number is the location of the peaks in Figs. 1 and 2 which 
correspond to the maximum of D\y and the shape of the functions of Fx around 
this peak arise as follows. The contributions to F\ from dY v + so ft and dT hard have 
opposite signs. The virtual and soft corrections start growing steeply at a slightly 
lower value of xi than the term from "hard" gluon radiation and the former level 
off slightly earlier, at Xi ~ 0.218, where the increase of dYhard is still rather steep. 
We have checked numerically that this does not depend on the choice of the (small) 
separation parameter x m i n . 

In the fits we use the generic form 

A 8 
Fx = 6(0.218 - x) {x _ 021 ° 8) 2 + Ai + ©(* - 0-218) £ B k x k , (9) 

where Q(x) denotes the step function. The coefficients A k ,Bk that specify the re- 
spective functions are given in Table 1. The corresponding functions for the neutrino 
distribution are obtained with the help of (JHJ). 

The polarisation dependent part of the neutrino distribution, F us = Fq s + 
Cf^-Fi s , is the term in the above distributions that is most sensitive to a cou- 
pling k r [T3^ . especially for neutrino energies x v < 0.6. For instance, for kr = 0.1 
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the function F us deviates from its SM value by about 4% for 0.25 < x v < 0.6, and 
for x v < 0.25 the deviations are even larger. In order to obtain maximal sensitivity 
to k,r in future data analyses one may use likelihood functions that can be derived 
[HI] from © and dU). 

Finally we consider the possibility that Im(a*/3) ^ 0. This happens in a natural 
way in SU{2) L x SU{2) R x U{1) extensions of the SM, where the Higgs sector is 
sufficiently complicated to allow for a CP-violating phase in the L — R gauge boson 
mixing matrix already at tree level. After diagonalisation of this mass matrix this 
phase is transported to the charged current interaction (j2J). An observable to check 
for such a phase is the expectation value of the T-odd triple product O = pV (p& x s), 
where are the directions of flight of the b quark and the lepton, respectively, in 
the top quark rest frame. (At this order in the QCD coupling the decay amplitude 
has no absorptive part that could generate < O >^ 0.) We have to order a s \ 

<0>= ^-( J OdT B+v+soft + J OalT hard ), (10) 

where T SL denotes the order a s partial decay width for (JTJ). Because the deviations 
of the moduli of the parameters a and (3 from their SM values are not expected to be 
large we use T^ 1 for the normalisation in (jlUJI . We use Pf^f = 0.168 GeV and 0.153 
GeV to lowest order and to order a s , respectively. We get < O >= clm(a*/3)S, 
where c = 0.0050 at tree level and c = 0.0042 including the order a s corrections. 
This implies that the effect is too small in order to obtain an interesting sensitivity 
to Im(a*/3) (i.e. at the percent level) from semileptonic top quark decays. 

In summary we have computed, for a tbW vertex with left- and right-handed com- 
ponents, double differential lepton distributions to order a s for polarised semilep- 
tonic top decay, and we presented the QCD corrections in terms of compact parame- 
terisations. These formulae should be a useful module in the theoretical description 
of t and/or t production and decay at next-to-leading order in a s , including non- 
standard top couplings. 
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FIGURE CAPTIONS 

Fig. 1. The order a s contributions to the polarisation independent part of the 
lepton distributions (jSJ), ©• 

Fig. 2. The order ot s contributions to the polarisation dependent part of the lepton 
distributions ©, ©• 
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Figure 1: 
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